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BESCHREDBUNG 

Rekonstruktion des Stromungsflusses in einem GefaBsystem 

Die Erfindung betrifit eine Datenverarbdtimgseinrichtung und ein Verfahien zur Rekon- 
struktion des Stromungsflusses in einem GefaBsystem sowie eine Anordnung, die eine derar- 
5 tige Datenverarbdtungseinrichtung enthalt 

Aus der US 2002/0041654 Al ist eine Rontgeneinrichtung bekannt, mit welcher die Aus- 
bieitung einer Kontrastmittelinjektion im GefaBsystem eines Patienten beobachtet werden 
kann. Der Anwender muss dabei die Rontgenapparatur wahiend der Kontrastmittekus- 
10 breitung bewegen, urn die Ausbreitungsfront im Bild zu behaltea Um diesen Bewegungs- 

vorgang ohne Strahlenbelastung fur den Patienten uben zu konnen, wird mit HDfe gespeicher- 
ter Aufriahmen eine simulierte Kontrastmittelausbreitung auf einem Monitor daigestellt Eine 
ModeUierung des realen Stromungsflusses ist bei dem System nicht vorgesehen. 

15 Vor diesem EHntergrund war es eine Aufgabe der vorliegenden Erfindung, Mittel zur Rekon- 
struktion des Stromungsflusses in einem GefaBsystem bereitzustellen, welche eine moglichst 
prazise ModeUierung der realen Verhaltnisse ermoglichen. 

Diese Aufgabe wird durch eine Datenverarbdtiingseinrichtung mit den Merkmalen des An- 
20 spruchs 1, durch eine Anordnung mit den Merkmalen des Anspruchs 8 sowie durch ein Ver- 
fahien mit den Merkmalen des Anspruchs 10 gelost Vorteilhafte Ausgestaltungen sind in den 
Unteranspriichen enthalten. 

Die erfindungsgemaBe Datenvea^beitxmgseinrichtung dient der Rekonstiuktion des Stro- 
25 mungsflusses in einem GefaBsystem. Das GefaBsystem kann insbesondeie ein interessierender 
Teil der BlutgefaBe sein, z.B. das System der HerzkranzgefaBe, ohne dass die Erfindung je- 
doch auf diesen Anwendungsfall beschrankt ware. Die Datenveiatbeitungseinrichtung umfasst 
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einen Speicher, welcher Messdaten enthalt, die eine beobachtete fortschreitende Ausbreitung 
eines (ersten) Mediums im genannten GefaBsystem beschreiben. Bei dem Medium kann es 
sich prinzipiell umjedes Objekt handeln, das eine fortschreitende Bewegung durch das GefaB- 
system ausfuhrt Insbesondere kann das Medium ein im GefaBsystem stromendes Fluid wie 
5 zum Beispiel eine Kontrastmittelinjektion sein. Die Beschreibung der Ausbreitung des 
Mediums kann durch den gemessenen zeitlichen Verlauf der Position des Mediums, zum 
Beispiel der Position einer Kontrastmittelfront oder eines Kontxastmittel-Konzentrations- 
maximums, erfolgen. Die Datenverarbeitungseinrichtung ist dazu eingerichtet, aus den 
genannten Messdaten eine modellierte Ausbreitung beziehungsweise "Modell-Ausbreitung" 
10 eines zweiten Mediums im GefaBsystem zu rekonstruienen. Das zweite Medium kann mit dem 
ersten identisch oder hiervon verschieden (z.B. Kontrastmittel vs. Blut) sein. Bei der Rekon- 
struktion der 'Modell-Ausbreitung" sollen fur das betrachtete GefaBsystem die folgenden zwei 
Bedingungen erfullt werden: 



15 - Erstens soli die Differenz zwischen der beobachteten Ausbreitung und der Modell- 
Ausbreitung minimal sein. Mit anderen Worten soil sich die gesuchte Modell- 
Ausbreitung moglichst eng an die Messdaten anlehnen. 



Zweitens soli die ModeU- Ausbreitung in Bezug auf das GefaBsystem monoton 
20 f ortschreitend sein. Das bedeutet, dass sich ein Medium entsprechend der Modell- 
Ausbreitung nur "vorwarts" in eine Richtung durch das GefaBsystem bewegen kann. 



Mit der beschriebenen Datenverarbeitungseinrichtung ist es moglich, aus Messdaten iiber die 
Ausbreitung eines (ersten) Mediums in einem GefaBsystem konsistent und in stabiler Weise ein 
25 Modell zu rekonstruieien. Dies gelingt dadurch, dass eine Modellfunktion global, das heiBt das 
gesamte betrachtete GefaBsystem beriicksichtigend an die Messdaten angepasst wird, wobei 
gleichzeitig die notwendige Randbedingung der monotonen Ausbreitung beachtet word Wenn 
beispielsweise die Ausbreitung einer Kontrastmittelinjektion im GefaBsystem eines Patienten 
rnit Rontgenaufhahmen gemessen wird, ergeben sich in der Praxis aufgrund von zahlreichen 
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StQreinflUssen stark streuende Daten. Dabei kann es insbesondere sein, dass keine monotone, 
das heiBt stetig fortschredtende Ausbreitung der Injektion durch die Daten beschrieben wird, 
sondem scheinbar Spriinge und/oder Ruckflusse auftreten. Ln Gegensatz hierzu wird bei der 
Berechnung der Modell- Ausbreitung fur das gesamte GefaBsystem die Monotonie strong 
5 beachtet, wodurch einerseits eine zwingende Randbedingung erfullt und andererseits eine 
signifikante Stabilisierung der Stromungsrekonstruktion erreicht wird 

GemaB einer ersten Weiterbildung der Datenverarbeitungseinrichtung fuhrt diese die 
Rekonstruktion derart aus, dass die Modell- Ausbreitung zusatzlich einen moglichst glatten 
10 Verlauf hat Wird die Modell- Ausbreitung allgemein als eine Funktion f(T) einer Variablen x 
beschrieben, so kann der "glatte Verlauf 1 insbesondere dadurch definiert werden, dass der 
Betrag der zweiten Ableitung, If I, im Mittel minimal sein soil Wenn f(x) beispielsweise die 
Position einer Ausbreitungsfront zur Zeit % ist, dann entspricht die Glattheit anschaulich einer 
moglichst beschleunigungsarmen Ausbreitung, 

15 

GemaB einer weiteren bevorzugten AusMmmgsform der Daterrverarbeitmgseinrichtung 
enthalt der Speicher als Messdalen sogenannte Bolusankunftszeiten m. Dabei ist das 
GefaBsystem logisch in einzelne Abschnitte unterteilt, die durch die Indices i=l,.JNf gekenn- 
zeichnet sind, und die Bolusankunftszeit m ist definiert als die in einer Messung ermittelte Zeit, 

20 welche ein (erstes) Medium ausgehend von einem vorgegebenen Startpunkt fur das Erreichen 
des GefaBabschnittes mit dem Index i benotigt Die dieser Ausgestaltungsform zugrunde 
liegende Diskretisierung des Problems ist einerseits fur die Behandlung in einer digitalen Da- 
tenverarbeitung^einrichtung gut geeignet, andeterseits entspricht sie auch der Tatsache, dass 
die beobachtete Ausbreitung eines Mediums in der Regel durch Messungen zu diskreten Zeit- 

25 punkten abgetastet wird. 



In Weiterbildung der vorstehend beschriebenen Ausfuhrungsform ist die Datenverarbdtongs- 
einrichtung femer dazu eingerichtet, Modell-Bolusankunftszeiten fur die GefaBabschnitte i so 
zu berechnen, dass die folgenden beiden Bedingungen erfullt sind: 
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(1) 



und 



£ = £|m r - istminimaL 



(2a) 



5 Gledchung (1) bringt die Forderung nach einer monoton fartschreitenden Ausbreitung zum 
Ausdruck, wahiend (2a) gjobal fur das GefaBsystem tine moglichst enge Anlehnung der 
ModeU-Bolusankujoftszeiten ti an die beobachteten Bolusanknnftszeiten m, fordert Dabei 
geben die Werte p(i) jeweils den Index des GefaBabschnittes wieder, welcher in Stromungs- 
richtung gesehen vor dem GefaBabschnitt mit dem Index i liegL 



Weiterhin kann die vorstehend beschriebene Datenverarbeitimgseiiirichtung derart erganzt 
werden, dass zusatzlich die Forderung 



beriicksichtigt wird, wobei die Menge I die Indices aller "inneren" GefaBabschnitte rrdteinem 



15 Vorganger und einem Nachfolger enthalt und t { die diskrete Approximation der zweiten 
Ableitung im GefaBabschnitt i ist Durch Bedingung (2b) wird die Forderung nach einem 
moglichst glatten Vedauf der ModeU- Bolus ankuiiflszeiten U ausgedriickt Bei der Kombination 
der Forderungen (2a) und (2b) wild in der Regel ein Gewichtungsfaktor anzusetzen sein, 
durch dessen Festlegung die relative Bedeutung der beiden Kriterien nach Wahl gewichtet 

20 werden kann. 

Bei einer Weiterbildung der vorstehend definierten Ausfuhrungsfoim ist die Datenverarbei- 
tungseinrichtung dazu eingerichtet, die Modell-Bolus ankunftszeiten ti mit Hilfe des Algorithmus 
der linearen Programmierung zu berechnen. In der Fonnulierung gemaB den Gleichungen (1), 
25 (2a) und gegebenenfalls (2b) lasst sich die gesuchte ModeU-Bolusankunftszedt namlich als 
Losung eines linearen Optiinierungsproblems darstellen, 



10 




ist minimal 



(2b) 
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so dass auf hierfur entwickelte Algorithmen der linearen Programmierung wie z3. das 
Simplexverfahren zuruckgegriffen werden kann. 

5 Vom GefaBsystem, in welchem der Stromungsfluss rekonstruiert werden soli, kann im 
Speicher insbesondere ein dieidimensionales Modell abgelegt sein. Dreidimensionale 
GefaBmodelle sind im Rahmen medizirrischer Anwendungen zunehmend anzutreffen und von 
grofiem Nutzen bei der Dateninterpretation durch den behandelnden Arzt. 

10 GemaB einer anderen Weiterbildung der Datenverarbdlxmgseinrichtxing kann diese dazu 
eingerichtet sein, bei der Rekonstruktion der Modell- Ausbreitung physiologische Randbe- 
dingungen zu beriicksicMgen. Beispielsweise kann bei einer Kontrastmittelinjektion die 
Erhaltung der injizierten Kontrastmittelmasse bei deren Ausbreitung im GefaBsystem beriick- 
sichtigt werden, wenn die zur Beschreibung der Ausbreitung verwendeten Messdaten be- 

15 ziehungsweise die gesuchte Modell- Ausbreitung eine Beriicksichtigung des Massenwertes 
vorsehen. Durch die Berucksichtigung physiologischer Randbedingungen, die zwangslaufig 
erfullt sein miissen, kann die berechnete Modell- Ausbreitung gegeniiber StSreinfliissen weiter 
stabilisiert werden. 

20 Die Datenverarbeiteuagseinrichtung ist femer vorzugsweise mit einer Anzeigeeinrichtung wie 
etwa einem Monitor gekoppelt, urn doit die berechnete Modell- Ausbreitung graphisch dar- 
stellen zu konnen. Auf diese Weise konnen die Exgebnisse der Berechnung zum Beispiel einem 
behandelnden Aizt bei einer Katheteruntersuchung visuell zur Vetfiigung warden. 

25 Die Erfindung betriffi femer eine Anordnung zur Beobachtung des Stromungsflusses in einem 
Gefafisystem, wobei die Anordnung eine bildgebende Einrichtung zur Erzeugung von Abbil- 
dungen des Gefafisystems enthalt, und wobei aus den genannten Abbildungen Messdaten 
gewonnen werden konnen, die die fortschreitende Ausbreitung eines Mediums im GefaB- 
system beschreiben. Femer enthalt die Anordnung eine Datenverarbeitungs einrichtung der 
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oben erlauterten Art zur Rekonstruktion des Stromungsflusses in dem GefaBsystemu Das 
heiBt, dass die Datmveraibeitungseinrich^ einen Speicher mit Messdaten umfasst, die die 
beobachtete f ortschreitende Ausbreitung eines (ersten) Mediums beschreiben und die durch 
die bildgebende Einrichtung gewonnen wurden. Weitethin ist die Dateiwerarbedtungsein- 
5 richtung dazu eingerichtet, global fur das betrachtete GefaBsystem eine Modell-Ausbreitung 
des (zweiten) Mediums zu beiechnen, die die Differenz zwischen der beobachteten und der 
modellierten Ausbreitung minimiert und gjeichzeitig eine monotone Ausbreitung beschreibt 



Bei der bildgebenden Einrichtung kann es sich insbesondere urn einen Rontgenapparat han- 
10 deln, vorzugsweise urn eine rotierende Rontgeneinrichtung zur Erzeugung dreidimensionaler 
Rontgenaufhahmen ednes Korpers. In der Regel stehen derartige Gerate in Katheterlabors 
ohnehin bereit und konnen daher problemlos fur die Ermitdung der Ausbreitung einer Kon- 
ttastmittelinjektion und damit der Stromungsverhaltnisse im GefaBsystem eines Patienten 
eingesetzt werden. 

15 

Die Erfindung betrifft femer ein Verfahien zur Rekonstruktion des Stromungsflusses in einem 
GefaBsystem, welches die folgenden Schritte umfasst: 

a) Die Gewinnung von Messdaten, die eine beobachtete fortschreitende Ausbreitung 
eines Mediums im GefaBsystem beschreiben. 
20 b) Die Rekonstruktion einer ModeU- Ausbreitung ednes Mediums in dem GefaBsystem 
derart, dass 

die Differenz zwischen der beobachteten Ausbreitung und der Modell-Ausbreitung 
minimal ist, und 

die Modell-Ausbreitung monoton fortschreitend ist 

25 

Das Verfahien implementiert in allgemeiner Form die mit der oben beschriebenen Daten- 
verarbtitungseinrichtung beziehungsweise Anordnung ausfuhrbaren Schritte. Zur Erlauterung 
des Verfahrens und der hierzu moglichen Varianten sei daher auf die obige Beschieibung 
verwiesen. 
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Ln Folgenden wild die Erfindung mit Hilfe der beigefiigten Figuren beispielhaft erlautert Es 
zeigt 

Fig. 1 schematisch die Komponenten einer erfindungsgemaBen Rontgenanordnung zur 
5 Rekonstruktion des Stromungsflusses in einem GefaBsystem; 

Fig. 2 die Definition von Variablen und Indices in Bezug auf das dieidimensionale Modell 

eines GefaBbaumes; 
Fig. 3 einige der fur die lineaie Programcaierung benotigten Gleichungen. 

10 In Rgur 1 ist als stellvertretendes Beispiel fur den Einsatz der vorliegenden Erfibnidung eine 

mediziriische Rontgeneinrichtung dargestellt Diese enthalt einen Rontgenapparat 1, welcher an 
einem C- Ami angebracht ist und zur Erzeugung von Rontgemufoahmen urn einen Patienten 2 
rotieren kann. Die aufgenommenen Daten werden zur Weiterverarbeitung an eine Datenver- 
afbdtungseinrichtung 3 weiteigeleitet. 

15 

Die Ronlgenapparatur 1 kann insbesondere dazu eingerichtet sein, dieidimensionale Rekon- 
struktionen eines KQrpervolumens zu eizeugen, urn z.B. die raumliche Struktur eines GefaB- 
systems zu untersuchen. Ferner ist versucht worden, zeitliche Ablaufe wie insbesondere die 
Stromung beziehungsweise die Ausbreitung eines Fluids im GefaBsystem zu rekonstruieren 

20 (vgL H. Schmitt, M. Grass, V. Rasche, O. Schramm, S. Haehnel, K. Sator: "An x-ray-based 
method for the determination of the contrast agent propagation in 3-D vessel structures", IEEE 
Trans. Med. Imag., 21(3):251-262, 2002). Derartige "4D-Abbildungen" begegnen erheb- 
lichen technischen Problemen, beispielsweise der Notwendigkeit hoher Bildralen zur Er- 
fassung des Stromungsflusses, der Stabilitat der Volumensegmentierung (das heiBt Zuordnung 

25 von Voxeln zu Objekten) und der Bewegungskorrektur. Dariiber hinaus stellt sich das grund- 
satziiche Problem einer guten Stromungsrekonstruktion fur ein gesamtes GefaBsystem. Beim 
Stand der Technik werden diesbezuglich aus den Rontgenaufioahmen fur den GefaBbaum 
sogenannte Bolvisankunftszeiten OBAT: bolus arrival times) bestimmt und dem GefaBbaum 
zugsordnet Die dabei erhaltenen Daten zeigen jedoch eine starke Stteuung, welche unter 
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andeiem auf Rauschen in den verwendeten Rontgenaufiiahmen, auf Ungenauigkeiten bei der 
Segmentierung, auf einer unvollstandige Bewegungskorrektur und auf einer nicht exakten 
Bestimmung der Bolus ankunftszeiten beruhen. 

5 Vor diesem Hintergrund soil nachfolgend unter Bezugnahme auf die Figuren 1 bis 3 ein Ver- 
fahren fur eine stabile und zuverlassige Rekonstruktion der Stromungsvethaltnisse in einem 
GefaBsystem erlautert warden. Das Verfahien beginnt mit einer dxeidimensionalen Rekon- 
struktion des GefaBbaumes 6 basierend auf rotatorischen 3EX-R6ntgenaufliahmen. Fwelheiten 
dieser Rekonstruktion konnen der einschlagigen Literatur entnommen werden (beispielsweise 

10 M. Grass, R. Koppe, E. Klotz, R. Proksa, M. H. Kuhn, H. Aerts, J. Op de Beek, R. Kem- 
kers: "3D reconstruction of high contrast objects using C-atm image intensifier projection 
data", Comput Med Imaging Graph., 23:311-321, 1999; M. Grass, R. Koppe, R Klotz, J. 
op de Beek, R. Kemkers: "3D reconstruction and imaging based on C-aim systems", Med. 
BioL Eng. Comp., 37(2): 1520- 1521, 1999). Nach der Segmentierung der GefaBe ist die 

15 Struktur des GefaBbaumes 6 bekannt 

Der GefaBbaum 6 kann gemaB Figur 2 in einzelne. Abschnitte oder Cluster aus Volumen- 
elementen (Voxeln) unterteilt werden, welche durch einen Index i=l,...N bezeichnet werden. 
Per definitionem sei dabei der groBte Index N dem Startpunkt des GefaBbaumes zugeordnet 
20 Des Weiteren gebe die Funktion p(i) jeweils den Index des GefaBabschnittes wieder, welcher 
in Stromungsrichtung gesehen unmittelbar vor dem GefaBabschnitt mit dem Index i liegt Diese 
Funktion wird insbesondere bei Verzweigungspunkten des GefaBbaumes benotigt, da dort 
keine in Enerschritten forflaufende Nummerierung beider GefaBzweige moglich ist 

25 Nachdem der GefaBbaum mit Hilfe einer ersten Serie von angiografischen Rontgenaufioahmen 
eimittelt worden ist, erfolgt eine zweite Aufhahmeserie, bei welcher die Rontgenapparatur 1 
vorzugsweise feststeht, das heiBt nicht urn den Patienten 2 rotierL Wahrend der zweiten 
Aufhahmeserie wird - beispielsweise uber einen Katheter - ein Kontrastmittel in das GefaB- 
system des Patienten 2 injiziert, und die resultieiende Ausbreitung des Kontrastmittels im 
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GefaBsystem wird auf den erzeugten R6ntgerpx)jektionen verfolgt Mit Hilfe von Methoden 
der Bildverarbeitung kann dabei aus den Projektionsaufhahmen fur jeden Abschnitt i des Ge- 
faBbaumes die zugehorige Bolus ankunftszeit rr\ ermittelt werderL Hierbei handelt es sich urn 
die Zeit, die das Kontrastmittel vom Zeitpunkt der Injektion bis zum Erreichen des GefaSab- 
5 schnittes i benotigt Die Bolus ankunftszeiten konnen der Beschreibung der Stroniungsver- 
haltnisse in den GefaBen dienen, und aus ihnen konnen zum Beispiel die lokalen Stromungs- 
geschwindigkeiten berechnet werden. Wie bereits erwahnt wurde, sind die gemessenen Bo- 
lusankunftszeiten n\ jedoch aufgrund zahlreicher Storeinfltisse einem hohen Rauschen unter- 
worfen, Haufig verietzen sie daher grundlegende physiologische Randbedingungen wie 
10 beispielsweise die Monotonie der fortschreitenden Stromungsausbreitxmg. Das heiBt, dass ein 
in Stromungsrichtung weiter hinten gelegener GefaBabschnitt eine kurzeie gemessene Bolus- 
ankunftszdt m haben kann als ein vor ihm gelegener Abschnitt 

Zur Behebung dieser Problematik und zur Erzeugung einer stabilen Rekonstruktion der Str5- 
15 mungsverhaltnisse im GefaBsystem wird die Berechnung einer Modell- Ausbreitung in Form 
von ModeU-Bolus ankunftszeiten U vorgeschlagen, die im Rahmen einer globalen Optimierung 
die Monotoniebedingung erfullen und gleichzedtig moglichst dicht bei den Messdaten liegen. 
Mit der oben eingefuhrten NomenMatur kann die Monotoniebedingung fur die gesuchten 
ModeU-Bolus ankunftszeiten ti duich die folgende Gledchung ausgedriickt werden: 



20 



(1) 



Unter Beachtung dieser Randbedingung soil nun die Kostenfunktion 



25 




(2a) 



minimiert werden, welche die iiber den gesamten GefaBverlauf berechnete Abweichung 
zwischen den gemessenen und den modellierten Bolus ankunftszeiten quantifizierL Die 
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Berechnung der gesuchten GroBen ti kann dabei mit HBlfe der linearen Progranunierung 
erfolgen (vgl. D. G Lunenbergen 'linear and Nonlinear Progra mm i ng ", Addison- Wesley, 
Reading, 2. Auflage, 1984; W. H. Press, S. A, Teukolsky, W. T. Vetterling, B. P. Hannery: 
"Numerical Recipes in C'\ Kapitel 10.8, University Piess, Cambridge, 2. AufL, 1992). Bei 
5 der sogenannten Normalform dieser Technik geht es darum, die Funktion 

F = c-a (3) 
unter Beachtung der Randbedingungen 

S a=b (4) 
a t >0 (5) 

10 zu minimieren, wobed b , c und S gegebene Vektoren beziehungsweise eine gegebene Matrix 
sind und der Vektor a die gesuchte Losung darstellt 



Urn die Technik der linearen Programmierung fur die vorliegende Problemstellung einsetzen zu 
konnen, sollen zwei N-dimensionale Hilfsvektoren e + und £~~ sowie ein N-dimensionaler 

15 Hilftvektor A eingefuhrt und die folgenden Vektoren a,b 9 c sowie die Matrix S definiert 
werden: 



a = 



, A , 



(6) 



b = 



f m N 



(7) 



20 



c = 



'5. 



(8) 
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5 = 



^NxN ^NxN ^Nxlf OaTx(W-I) ^ 



(9) 



Dabei bezedchnen l M und 0 M M-dimensionale Vektoren, die eine Eins beziehungsweise eine 
Null in jeder Komponente enthalten; 1 M xm ist eine EMieitsmatrix der GroBe MxM; Omxl ist 
5 eine Matrix mit M Zeilen und L Spalten, deten Komponenten alle Null sind; und P ist die (N- 
l)xN Kopplungsmatrix, die definiert wird durch 



1 falls i = j 

-1 falls p(i) = j (10) 
0 sonst 



10 Da p(i) i , ist die Definition von P widerspnichsfirei. 

Unter Verwendung der vorstehend definierten GroBen soil durch lineare Programmierung die 
Funktion 

15 F = ?-a = £(e f + +er) (ID 



unter Beachtung der Randbedingungen 

f ( +e / + -e~ = m, fur i =1,..JV (12) 

20 t, - t m -A f =0 far i = U.N - 1 (13) 

f,:>0 (14) 

e, + £0 (15) 

ef;>0 (16) 

A,£0 (17) 
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minimiert werden. Gleichung (11) impliziert zusammen mit den Gleichungen (12), (14)- (17), 
dassfurjedenIhdexi=l v .,Nentweder e* oder e," Null sein muss. Ware namlich zum 
Beispiel > e* >0, so wifrdemitdemAnsatz e f ~ := e7 - 8*, e, + = 0 ein kleinerer Wert in 
5 Gleichung (1 1) erhalten. Da weiterhin 

8. = - e ~ = m i - t i (18) 
gilt, ist die Minimierung der GroBe F in Gleichung (11) aquivalent zur Minimierung von E in 
Gleichung (2a). In ahnlicher Weise ist die in Gleichung (1) definierte Neberibedingung aqui- 
valent zu den Gleichungen (13) und (17). Die mit Standardmefhoden der numerischen 
10 Mathematik ermittelbare Losung fur die Gleichungen (1 1)-(17) liefert somit zugleich die 
Losung des Minimierungsproblems von Gleichung (1) und (2a), falls die Messdaten hin- 
reichend grofi sind, urn die Neberibedingung U > 0 sicher zu steflen. Letzteres lasst sich jedoch 
in einf acher Weise durch Addition einer Konstanten zu alien Messdaten erreichen. 



15 Es sei darauf hingewiesen, dass es einen sogenannten zulassigen Basisvektor gibt, das heifit 
einen Vektor a , der alle Randbedingungen erfullt Mit der Definition 

x =max{m < \i =1,..JV} (19) 

kann der Vektor a namlich aus dan folgenden Komponenten zusammengesetzt werden 

F=(x,...<c) r (20) 
20 e + =(0 v ..0f (21) 

e - = (x - m p ...x - m N J (22) 
A = (0,...0) r . (23) 



Das oben erlauterte Verfahren kann zusatzlich im Rahmen einer sogenannten Regularisterung 
25 urn eine Kriterium erganzt werden, das einen glatten Verlauf der Losung fordert Nachfolgend 
soli diesbezuglich zusatzlich zu Gleichung (2a) die Summe der Absolutwerte der zweiten 
Ableitungpn der Bolusankunflszeiten ti minimiert werden. 
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Zu diesem Zweck wird ein weiterer Satz von Variablen 8 1 + ,..,8 n-m-x und 5 t ,...5 ^-m-i 
definiert Dabei steht M fur die Anzahl an GefaBabschnitten ohne Nachfolger, &h. der 
stromabwartigen Enden des betrachteten GefaBbaumes. Ohne Beschrankung der Allge- 
5 meinheit kann dabei eine Indizierung vorausgesetzt werden, bei der diese M GefaBabschnitte 
gerade die Indices N-M, ... N-l haben. Fur jeden Index i= 1, ... (N-M-l) soli S* -5 r die 
zweite Ableitung der Bolus ankunftszeiten beim GefaBabschnitt i xeprasentieren. Daher miissen 
diese Variablen mit dem Vektor der Bolus ankunftszeiten F durch die folgende Randbe- 
dingung vedcaiipft sedn: 



Dabei entspricht I Jl der Anzahl an GefaBabschnitten, die den GefaBabschnitt i als Vorganger 
haben. Abgesehen von Verzweigungspunkten des GefaBbaumes hat der GefaBabschnitt i nur 
einen Nachfolger j, so dass IJI = 1 ist und die rechte Seite von Gleichung (24) zur ttblichen 
diskreten Approximation der zweiten Ableitung kollabiert Falls die Lange der einzelnen 
15 GefaBabschnitte innedialb des GefaBbaumes erfieblich variiert, kann die Rinfuhrung von 
geometrischen Gewichtungen empfehlenswert sein. 

Nunmehr kann die lineare Programmierung dazu verwendet werden, die folgende GroBe zu 
minimietenz 



Ein kleiner Wert von Em bedeutet dabei, dass die Absolutwerte der zweiten Ableitungen der 
Bolusaokunftszeiten klein sind. Daher werden "glatte" Losungen t bevorzugt 



10 




(24) 



20 




(2b) 



Die Regularisierung erfbidert weiterhin die Einfuhrung eines Regularisierungsparameters X, 
25 welcher die beiden zu minimierenden Zielfunktionen, dh. den "Datentenn" E gemaB Glei- 
chung (2a) und den "Modellterm" I^n gemaB Gleichung (2b), gewichtet 
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Die fur die lineare Programmierung verwendeten GroBen konnen dann gemaB den Gleichun- 
gen (26) bis (29) von Figur 3 angesetzt werden, wobei zur Vereinfachung I=N-M- 1 ge- 
schrieben wird. Die zusatzJiche IxN Kopplungsmatiix C implementiert die rechte Seite der 
Bedingung von Gleichung (24) und wird gemaB Gleichung (30) von Figur 3 definieit 

5 

Bei der Regularisierung ist es in der Regel schwer, vorab einen "guten" Wert fiir den Regulari- 
sierungsparameter X zu wahlen. Zur Uberwindung dieses Problems wird vorliegend ein zwei- 
stufiges Verf ahren vorgeschlagen. Dabei wird in der ersten Stufe die Optdmierung ohne Regu- 
larisierung ausgefuhrt Dies liefert eine erste Schatzung F (0) fiir die Bolus anlainftszdten, mit 
10 denen sodann der "bestmogUche" Datenterm l£ 0) berechnet wird, dJu 

E io) =±fa-tr\ . (3D 

Nunmehr kann ein sinnvoller zu eneichender Wert fur E angesetzt werden, z.B. E= sE^ 0) mit 
s> 1. Da alles linear ist, kann die uispritagftche Mini^ 
15 E = s^ 0) ersetzt werden, was fur die lineare Programmierung zu den Gleichungen (32) bis (35) 
von Figur 3fuhrL 

GemaB Gleichung (36) ist es dabei weiterhin moglich, einen zulassigen Basisvektor zu finden. 
Gegebenenfalls ist es wiederum erforderlich, alle Werte ti (1) zu verschieben, falls das erhaltene 
20 t* (1) negaliv ist 

Altemadv kann im Ubrigen auch der "worst case" Modellterm basierend auf dem ersten 
Vektor F (0) berechnet werden und die Bedingung Em = sEm (0) mit s < 1 zum Mkrinrienings- 
problem des Datenterms E addiert werden. 

25 

SchlieBlich sei darauf hingewiesen, dass das vorgestellte Verfahren auch die Integration 
anderer physiologischer Ihformationen erlaubt, welche in Form von linearen Randbedingungen 
beschrieben werden konnen. Zum Bdspiel konnte die Masseneihaltung auf diese Weise for- 
muliert werden. 
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PATENTANSPRUCHE 



1. Datenverarbeitogseinrichtaig (3) zur Rekonstruktion des Stromungsflusses in einem 
GefaBsystem (6), enthaltend einen Speicher (4) mit Messdaten (mO, die eine beobachtete 
fortschieitende Ausbreitung ednes Mediums ini GefaBsystem (6) beschreiben, wobei die 
Datenverarbdtmgseimichtung (3) dazu eingerichtet ist, aus den Messdaten eine Modefl- 

5 Ausbreitung (U) eines Mediums im GefaBsystem derart zu rekonstruieren, dass fur das 
GefaBsystem 

die Differenz zwischen der beobachteten Ausbreitung und der Modell- Ausbreitung 
minimal ist, und 

die Modell- Ausbreitung monoton fortschieitend ist 

.0 

2. Datenverarbdtxingseinrichtung nach Anspruch 1, welche dazu eingerichtet ist, die Modell- 
Ausbreitung (ti) derart zu iekonstruieren, dass diese zusatzlich einen moglichst glatten Verlauf 
hat 



15 3. Datenverarbdtiingsdnrichtung nach Anspruch 1, 
dadurch ppkftriny eichnet 

dass der Speicher (4) als Messdaten Bolusankuriflszeiten rik enthalt, wobei die i= 1,...N 
Indizes fur verschiedene Abschnitte des GefaBsysterns (6) sind und eine Bolusankunftszeit mi 
die in einer Messung ermittelte Zeit ist, die ein Medium ausgehend von einem vorgegebenen 
20 Startpunkt fiir das Ikreichen des GefaBabschnittes i benotigt 
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4. Datenverarbdtimgseinrichtxing nach Ansprach.3, 
dadurch gekennzeichnet 

dass sie dazu eingerichtet ist ModeU-Bolusankunftszeiten ft) fur die GefaBabschnitte i derart 
5 zu berechnen, dass 

a,=*,-»ko^ 0 vi=i,..j\r-i (D 

und die Kostenfunktion 

2? = £h-f,| (2a) 

minimal ist, wobei jeweils p(i) der Index des in Stromungsrichtung vor dem GefaBabschnitt i 
1 0 liegenden Gefafebschnittes ist. 

5. Datenvem:bdtmgseinrichtung nach Anspruch4, 
dadurch gekennzeichnet 

dass sie dazu eingerichtet ist, zusatzlich in der Kostenfunktion die GroBe 

15 



ial 



(2b) 



zu berucksichtigen, wobei I die Indices aller GefaBabschnitte mit einem Vorganger und einem 
Nachfolger enthalt und t t die diskrete Approximation der zweiten Ableitung im 
GefeBabschnitt i ist 

20 

6. Datenverarbeitungseinrichtung nach Anspruch4, 
dadurch gekennzeichnet 

dass sie dazu eingerichtet ist die Modefl-Bolusankunftszeit ft) mit Hilfe der linearen 
Progmmmierung zu berechnen. 

25 
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7. Datenverarbdtimgseinrichtung nach Anspmch 1, 
dadurch ^ ekennT eichnet 

dass sie mit einer Anzdgeeinrichtung (7) zur graphischen DarsteUung der ModeU- Ausbreitung 
5 gekoppelt ist 

8. Anordnung zur Beobachtung des Sfrbmungsflusses in einem GefaBsystem (6), enthaltend 
eine bildgebende Einrichtung (1) zur Erzeugung von Abbildungen des GefaBsystems (6), aus 
welchen die fortschreitende Ausbreitung ones Mediums beschreibende Messdaten (n*) 

10 gewonnen werden kbnnen, sowie eine Datenverarbdtungseinrichtung (3) nach mindestens 
einem der Anspruche 1 bis 7 zur Rekonstruktion des Stromungsflusses in dem GefaBsystem. 

9. Anordnung nach Anspruch 8, 
dadurch gekennzeichnet 

15 dass die bildgebende Einrichtung eine Rontgenappaiatur (1) ist 

10. Verfahren zur Rekonstruktion des Stromungsflusses in einem GefaBsystem (6), umfessend 
die Schritte 

a) Gewinnung von Messdaten (m), die eine beobachtete fortschreitende Ausbreitung 
20 eines Mediums im GefaBsystem (6) beschreiben, 

b) Rekonstruktion einer ModeU- Ausbreitung (tj) eines Mediums in dem GefaBsystem 
derart, dass 

die Etifferenz zwischen der beobachteten Ausbreitung und der ModeU- Ausbreitung 

minimal jst, und 

25 - ^eModeU-Ausbreitungmonotonfortschreitendist 
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ZUSAMMENFASSIJNG 

Rekonstruktion des Stromungsflusses in einem GefaBsystem 

Die Erfindung betrifft ein Verfahien zur Rekonstruktion des Stromungsflusses beziehungsweise 
der Bolusankunflszeiten in einem GefaBsystem (6). Fur die Abschnitte (i) des 
5 GefaBbaumes (6) weiden zum Beispiel im Rahmen einer Kontrastmittelinjektion Bolus- 
ankunflszeiten (mi) gemessea Basieiend auf diesen Messdaten werden dann mit EBife der 
lineaien Programmierung Modell-Bohisankunftszeiten ft) beiechnet, welche einerseits durch 
Minimierung der Funktion E = £|m f - 1\ moglichst dicht an den Messdaten liegen und 
andererseits durch Einhaltung der Randbedingung A, =t i -t li(l) >0 die Monotonie der 
10 Ausbreitung gewahrleisten, wobei p(i) der Index des GefaBabschnittes vor dem 

GefaBabschnitt i ist Vorzugsweise wird durch eine zusatzliche Minimierung von E m =X|''"| 
ein moglichst glatter Verlauf der Ausbreitung erzwungen. 



(Figur2) 

15 
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Fig. 2 
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else 
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(32) 



(33) 



c = 



On 
On 
On 
Sn-i 
ti 



(34) 
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Fig. 3 



